Experimental manifestation of the breakpoint region in the current-voltage 
characteristics of intrinsic Josephson junctions 
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The experimental evidence of the breakpoint on the current-voltage characteristics (IVCs) of the 
stacks of intrinsic Josephson junctions (IJJs) is presented. The influence of the capacitive coupling 
on the IVCs of E^S^CaC^Oy IJJs has been investigated. At 4.2 K, clear breakpoint region 
is observed on the branches in the IVCs. It is found that the hysteresis observed on the IVC 
is suppressed due to the coupling compared with that expected from the McCumber parameter. 
Measurements agree well with the results obtained by the theoretical model. 

PACS numbers: 74.50.+r, 74.72.Hs, 85.25.Cp, 74.81.Fa 



It is well known that layered Bi2Sr2CaCu20 y sin- 
gle crystals represent natural stacks of atomic intrinsic 
Josephson junctions with interlayer spacing of 1.55 nnifi^ 
Their physical properties are different from those of con- 
ventional Josephson junctions because there exists the 
inductive or capacitive coupling between adjacent junc- 
tions due to extremely thin superconducting Cu02 layer 
of 0.3 nm. 

The influence of the coupling effect on the IVCs has 
been studied both theoretically and experimentally. In 
the absence of an external magnetic filed, the capacitive- 
coupling may become important because the thickness of 
CuC"2 layers is comparable to the Debye charge screen- 
ing length t\d. From this point of view, the capasitively 
coupled Josephson junctions (CCJJ) model was proposed 
to describe the phase dynamics of IJJs^ and was used 
to analyze their IVC. On the other hand, it was also 
theoretically shown that the diffusion current plays an 
important role in the IJJ stack£ Therefore, the CCJJ 
model including the diffusion current (CCJ J+DC model) 
may be more suitable for the description of IJJs under 
no external magnetic field. Recently, Shukrinov et al. 
have been systematically studied the phase dynamics of 
IJJs by using the CCJJ+DC model and predicted that 
the capacitive-coupling effect induces the appearance of 
new fine structures, which is called as breakpoint region 
(BPR), together with an equidistant branch structure in 
the IVCiS^iS The BPR is caused by the parametric res- 
onance between Josephson and plasma oscillations, and 
is determined as a region between the breakpoint current 
Ib p , at which the resonance is happened and the longitu- 
dinal plasma wave (LPW) is excited, and the transition 
current Ij to another branch in the IVC. However, to our 
knowledge, the experimental observation of the BPR in 
the IVC of the IJJs has been not done yet. 

In this letter, we show the first experimental results 
which demonstrate the BPR in the IVC. 

Bi2Sr 2 CaCu20 y single crystals with critical tempera- 
tures T c of ~ 85 K were grown by a conventional melting 
method. The grown single crystals were glued on glass 



substrates and the cleaved in air to obtain their fresh sur- 
face. Subsequently, Au thin films of 20 — 40 nm thickness 
were deposited on the cleaved surfaces. And, the mesas 
with areas S of 9, 16, 25, 36, 49, 64 /zm 2 , consisting of 
5-20 IJJs, were fabricated by using electron beam lithog- 
raphy, photolithography and Ar ion milling. Their IVC 
along the c-axis direction were measured using a three 
terminal method. Here, we present the results for two 
samples: Nml-25 and Nm5-25 with S — 25 /j,m 2 . 

Figure QJa) shows the typical IVC of the fabricated 
mesa (Nml-25). The IVC exhibits a multiple branch 
structure with large hysteresis. From the number of 
branches, it is found that the mesa contains of 8 IJJs. 
The critical current I c of constituent IJJs is nearly equal 
to each other and is 640 /iA, corresponding to the cur- 
rent density of J c = 2160 A/cm 2 . The voltage spacing 
between branches at I = I c is ~ 30mV. FigureQJb) shows 
the enlarged view of the IVC in the low bias region. We 
note that the minimal currents of each branch, i.e., the 
value of current at jump point Ij, are not identical in con- 
trast to the good homogeneity of J c 's. This feature is con- 
sistent with the prediction of the CCJJ+DC model^^ 
The ununiformity of the I/s can also be caused by self- 
heating but the self-heating effect can be ignored for this 
sample because the voltage spacing between branches at 
I c is independent of the number of junctions in the volt- 
age state in the stack. Therefore, the observed ununi- 
formity is predominantly due to the coupling between 
adjacent IJJs in the stack. Furthermore, the branches 
demonstrate some structure before transition to another 
branch. A detailed structure before Ij, which was ob- 
tained by sweeping the bias current up and down repeat- 
edly, is shown in Fig. [ljc). This is very similar to the 
BPR predicted in Ref. 7. It also looks like the subgap 
structure caused by the resonance of the Josephson oscil- 
lations and phonon modes i 10 ' 11 ! 12 However, we can rule 
out this phonon excitation because it should appear at 
V ~ 3, 6, 8 mV and these voltages are inconsistent with 
our results. Therefore, we consider that the observed 
structure is the experimental evidence of the BPR. 
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FIG. 1: The IVC of the sample Nml-25 at 4.2 K. (a) Over- 
all IVC. (b) Same IVC in the low bias region, (c) IVC on 
expanded scales showing the BPR. 



As it was shown in Ref. 8, in the CCJJ+DC model the 
hysteresis of IVC depends on the coupling parameter a = 
eru/ (d s di) and the dissipation parameter /3 = 
where e is the dielectric constant, d s is the thickness of 
superconducting layer, dt is the thickness of insulating 
layer and j3 c = 2eCR 2 I c /h is the McCumber parameter 
(C : junction capacitance, R : junction resistance, h : 
Plank's constant). If a — (no coupling) or it is negli- 
gible small, the IVC is consistent with that obtained in 
the resistively and capacitively shunted junction (RCSJ) 
model. Therefore, we compare the hysteresis of the ex- 
perimental IVCs with results of the RCSJ model. How- 
ever, for IJJs at low temperature the simple RCSJ model 
with constant R might not be used because of strongly 
nonlinear IVC, as is seen in Fig. [IJa). Thus, the IVC 
at 77 K, which becomes almost linear in a small voltage 
range, was examined. 

Figure [Ha) shows the IVC of the mesa shown in Fig. 
[TJ measured at 77 K and Fig. [^b) shows its low bias 
current region. As we can see, the multiple branches 
are almost linear but the large hysteresis still remains. 
And it is found that even at 77 K 7 c 's of the IJJs in the 
stack are almost identical. In contrast, the Ij of each 
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FIG. 2: The IVC of the sample Nml-25 at 77 K. (a) Overall 
IVC. (b) IVC in the low bias region. 



branch is not uniform, as is seen in Fig. [3c). However, 
we couldn't observe the BPR at 77 K. This is maybe 
because the thermal fluctuation makes the resonance a 
delicate state. Therefore, in this CcLSG, I bp 

and Ij coincide 

with each other. In Fig. HJb), we indicate the Ibp on the 
outermost branch and the return current I r from the first 
branch. Because the outermost branch extrapolates to 
zero current at zero voltage, from the sum voltage jump 
of AV = 39.1 mV at I = I c {= 240 fiA) and N = 8 one 
can estimate that the junction resistance R of each IJJ 
is 20.4 £1 at 77 K. Additionally, the junction capacitance 
C — eS/d t is calculated to be 1.84 pF by assuming e r = 
10. Substituting these parameters, we obtain /3 C (77 K) w 
560 for the presented mesa. Therefore, for the RCSJ 
model (a — 0) we expect I^} CSJ ~ 13 /jA by using the 

relation I r /I c = 4/(7r/3~ 1/2 ) for f3 c » l. 13 This value is 
significantly different from the I r — 45 /iA obtained from 
Fig. [2](b). Such a difference between experimental I r 
and I? CSJ can be explained by the capacitive coupling 
as shown in Ref. 8. Although another possible reason for 
the difference is the influence of thermal fluctuations, 14 
it seems unlikely that the fluctuations lead to such large 
difference discussed above. 

Figure [3] shows the temperature dependence of I Cl hp, 
I r and Ibp — I r for a different mesa (Nm5-25) with S = 
25 fim 2 and N = 8. The I C (T) is slightly deviated from 
the Ambegaokar-Baratoff dependence due to the effect 
of thermal fluctuation on the switching current. The I r 
hardly depends on T at low temperature (T < 40 K) 
and gradually increases with increasing temperature up 
to 80 K. Such behavior is typical for IJJs i 15 ' 16 It is also 
found that the Ibp shows T dependence similar to the I r . 
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FIG. 3: Temperature dependence of I c , hp, Ir and h p — I r 
of the sample Nm5-25 {S = 25 fim 2 and N = 8). 
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FIG. 4: The simulated IVC of a stack of 10 IJJs at a = 1, /3 = 
0.2 and 7 = 0. The inset shows the enlarged BPR for the last 
two branches in IVC. 

Furthermore, we note that 1^ — I r is almost independent 
of T for T < 40 K but at higher temperatures it increases 
with T. As shown in Ref. 7, the Ib p depends on a, (3 



and the wave number of excited LPW. Assuming a is 
independent of T, the increase in h p — I r may arise from 
the change of the excited plasma mode. However, to 
prove this, a further experiment would be required. 

In the framework of the CCJJ+DC model, the phase 
dynamics of IJJs can be described by the system of equa- 
tions 

d 2 (fe T . r,d(pi , . 

—j— = I - sin tpi - £f— + a(sm <p i+1 + sin ip t _ x 

-2 B m< Pt )+ a p{— + — -2— j (1) 

for the gauge-invariant phase difference ifi(t) = 9e+i(t) — 
9t(t) — ^ J g +1 dzA z {z,t) between superconducting lay- 
ers in the stacks with TV IJJs. Here 6i is the phase of 
the order parameter in S-layer £, A z is the vector poten- 
tial in the barrier. Figure [4] shows the simulation result 
of the IVC for a stack of 10 IJJs at coupling parameter 
a = 1, dissipation parameter (3 — 0.2 and with nonperi- 
odic boundary conditions 7 = 0, where 7 = s/sq = s/sn 
and s, So, sjy are the thickness of the middle, first, and 
last superconducting layers, respectively. In this figure, 
the dashed line corresponds to the IVC for a = 0. It 
is found that Fig. 2] replicates the experimental obser- 
vations descried above, such as the equidistant branch 
structure, BPR, discrepancy of Ij and the suppression of 
the hysteresis. 

In summary, we have experimentally investigated the 
influence of the coupling between IJJs on the IVC of 
Bi2Sr2CaCu20 a . We observed clear breakpoint region on 
the branches in the IVC and found that the coupling leads 
to the reduction of the hysteresis of the IVC at higher 
temperatures. Our results are in a good agreement with 
recent theoretical work based on the CCJJ+DC model. 
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